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We present a search for a narrow resonance in the inclusive diphoton final state using ~ 2.7 fb~^ 
of data collected with the DO detector at the Fermilab Tevatron pp Collider. We observe good 
agreement between the data and the background prediction, and set the first 95% C.L. upper limits 
on the production cross section times the branching ratio for decay into a pair of photons for 
resonance masses between 100 and 150 GeV. This search is also interpreted in the context of several 
models of electroweak symmetry breaking with a Higgs boson decaying into two photons. 

PACS numbers: 14.80.Bn, 13.85.Rm, 13.85.Qk 



At a hadron collider, diphoton (77) production allows detailed studies of the Standard Model (SM) as well 
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as searches for new phenomena, such as new heavy res- 
onances 0, extra spatial dimensions or cascade decays 
of heavy new particles Within the SM, continuum 
77 +X production is characterized by a steeply- falling 77 
mass {Mjy) spectrum, on top of which a heavy resonance 
decaying into 77 can potentially be observed. In particu- 
lar, this is considered one of the most promising discovery 
channels for a light SM Higgs boson at the LHC ;1, de- 
spite the small branching ratio of B{H 77) ^ 0.2% for 
110 < Mh < 140 GeV At the Tevatron, the dom- 

inant SM Higgs boson production mechanism is gluon 
fusion {gg H, or GF), followed by associated produc- 
tion with &W or Z boson [qq' VH, V = W,Z), and 
vector boson fusion [VV H, or VBF) 0,i,[9|. While 
the SM Higgs production rate at the Tevatron is not suf- 
ficient to observe it in the 77 mode, the Hgg and H^^ 
couplings, being loop- mediated, are particularly sensitive 
to new physics effects. Furthermore, in some models be- 
yond the SM % B{H 77) can be enhanced signifi- 
cantly relative to the SM prediction. 

In this Letter, we present a search for a narrow reso- 
nance in the Mj^^ spectrum using a data sample collected 
by the DO detector [l3| at the Fermilab Tevatron Col- 
lider. The selection of an inclusive jj+X sample and the 
use of the M^^ spectrum make the results of this search 
quasi- model independent. We use the SM Higgs boson 
(_ffsM) with H — !■ 77 as a reference model, resulting in 
the first such search at the Tevatron, and a forerunner 
to similar planned searches at the LHC. Additionally, we 
consider other models of electroweak symmetry breaking 
(EWSB) with enhanced B{H 77) as a consequence 
of suppressed Higgs couplings to either (i) all fermions 
(known as bosonic or fermiophobic Higgs boson, H{); (ii) 
only down- type fermions (^^u, as it gives mass to up- type 
fermions); or (iii) only top and bottom quarks (known as 
electroweak Higgs boson, i?ew)- AH models considered 
have SM-like production cross sections, with the excep- 
tion of the Hf and i/ew models, where GF is absent or 
has a negligibly small cross section. 

The subdetectors most relevant to this analysis are: 
the central tracking system, composed of a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT) 
embedded in a 2 T solenoidal magnetic field, the cen- 
tral preshower (CPS), and the liquid-argon and uranium 
sampling calorimeter. The CPS is located immediately 
before the inner layer of the calorimeter and is formed 
of one radiation length of absorber followed by several 
layers of scintillating strips. The calorimeter consists of 
a central section with coverage of I77I < 1.1 [lH, and two 
end calorimeters covering up to \r]\ ~ 4.2. The electro- 
magnetic (EM) section of the calorimeter is segmented 
into four longitudinal layers (EMz, i — 1,4) with trans- 
verse segmentation of A77 x = 0.1 x 0.1 ll|, except 
in EMS, where it is 0.05 x 0.05. The calorimeter is well- 
suited for a precise measurement of electron and photon 
energies, providing a resolution of ~ 3.6% at energies of 



'-^ 50 GeV. The data used in this analysis were collected 
using triggers requiring at least two clusters of energy 
in the EM calorimeter and correspond to an integrated 



luminosity of 2.7 ± 0.2 fb'^ [12|. 



Events are selected by requiring at least two photon 
candidates with transverse momentum px > 25 GeV 
and \ri\ < 1.1, for which the trigger requirements are 
fully efficient. The photon pT is computed with re- 
spect to the reconstructed event primary vertex (PV) 
with the highest number of associated tracks, which is 
required to be within 60 cm of the geometrical center 
of the detector along the beam axis. The PV recon- 
struction efficiency in + X events is ^ 98%, with 
~ 95% probability to match the true vertex. Photons 
are selected from EM clusters reconstructed within a 
cone with radius TZ = y/ (A?])^ -I- (Ac/))^ — 0.2 by re- 
quiring: (i) > 97% of the cluster energy is deposited in 
the EM calorimeter; (ii) the calorimeter isolation vari- 
able / = [£;tot(0.4) - £;em(0.2)]/Sem(0.2) < 0.1, where 
-Etot('7^) (^^em(^)) is the total (EM) energy in a cone of 
radius TZ; (iii) the energy-weighted shower width in the 
r — plane in EM3 is < -^14 cm; and (iv) the scalar sum 
of the PT of all tracks {p^^^J originating from the pri- 
mary vertex in an annulus of 0.05 < TZ < 0.4 around the 
cluster is < 2 GeV. To suppress electrons misidentified as 
photons, the EM clusters are required to not be spatially 
matched to tracker activity, either a reconstructed track, 
or a density of hits in the SMT and CFT consistent with 
that of an electron [l3|. To suppress jets misidentified as 
photons, a neural network (NN) is trained using a set of 
variables sensitive to differences between photons and jets 
in the tracker activity and in the energy deposits in the 
calorimeter and CPS: p^^]^, the numbers of cells above 
a threshold in EMI within 7^ < 0.2 and 0.2 < 7^ < 0.4 
of the EM cluster, the number of CPS clusters within 
TZ < 0.1 of the EM cluster, and the squared-energy- 
weighted width of the energy deposit in the CPS. The 
NN is trained using 77 and dijet Monte Carlo (MC) sam- 
ples and its performance is verified using a data sample 
of Z — > £+£"7 {£ — e, /i) events. Figure [T^ compares the 
NN output (Onn) spectrum for photons and jets. Pho- 
ton candidates are required to have Onn > 0.1, which 
is ~ 98% efficient for real photons and rejects ~ 50% of 
misidentified jets. Finally, M^^, computed from the two 
highest PT photons, is required to be > 60 GeV. In total, 
5608 events are selected in data. 

All MC samples used in this a naly sis are generated 
using PYTHIA [ill with CTEQ6L [l^ parton distribu- 
tion functions (PDFs), and processed through a GEANT- 
based [l3| simulation of the DO detector and the same 
reconstruction software as the data. Signal samples 
are generated separately for GF, VH and VBF pro- 
duction and normalized using the theoretical cross sec- 
tions 0, H, 0] and branching ratio predictions from hde- 
CAY [g. 

This analysis is affected by instrumental backgrounds 
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FIG. 1: (color online), (a) Normalized Ojvjv spectrum for photons and jets, (b) M~f^ spectrum in data (points) compared to 
the total background parameterization (solid line), including the DDP contribution derived via a sideband fit, and the total 
background prediction (dashed line) inside the search region for Mh = 130 GeV. The inset figure compares the data to the 
total background prediction inside the search region including its one standard deviation (s.d.) uncertainty band, as well as 
the expected Hqm signal scaled by a factor of 50. (c) Observed and expected 95% C.L. upper limits on cr x _B as a function of 
Mh- Also shown are the predictions for ct x _B in the different EWSB scenarios discussed in the text. 



such as 7+jet, dijet and Z/7* e+e^ (ZDY) produc- 
tion, with jets or electrons misidentified as photons, as 
well as an irreducible background from direct diphoton 
production (DDP). All backgrounds, except for ZDY, are 
estimated directly from data. 

The ZDY background is estimated using the MC sim- 
ulation, normalized to the next-to-next-to-leading-order 
cross section The selection efficiencies determined 

by the MC simulation are corrected to the corresponding 
values measured in the data. On average each electron 
has a 2% probability to satisfy the photon selection cri- 
teria, mainly due to the inefficiency of the track-match 
veto requirements. The total contribution from ZDY is 
estimated to be 88 ± 10 events. 

Backgrounds due to 7-l-jet and dijet events are di- 
rectly estimated from data by using a 4 x 4 matrix back- 
ground estimation method [l8|. After final event selec- 
tion, a tightened Onn requirement {Onn > 0.75) is 
used to classify the events into four categories depending 
on whether the two highest-pr photons, only the lead- 
ing photon, only the trailing photon or neither of the 
two photons, satisfy this requirement. The correspond- 
ing numbers of events, after subtraction of the estimated 
ZDY contributions, are denoted as Npp, Npj, Njp and 
Nff. The different relative efficiency of the O^n > 0.75 
requirement between real photons and jets allows the 
estimation of the sample composition by solving a lin- 
ear system of equations: {Npp, Npf , Nfp, Nf f)'^ = £ x 
{Nj^jN^j, Njj, Njj)'^ , where N^^ {Njj) is the number of 
77 (dijet) events and N-yj (Njy) is the number of 7+jct 
events with the leading (trailing) cluster as the photon. 
The 4x4 matrix £ contains the efficiency terms (pa- 
rameterized as a function of \ri\), estimated in photon 
and jet MC samples and validated in data. The esti- 
mated sample composition is Nj~^ = 3155 ± 125 (stat), 
N-^j+j-( = 1680 ± 149 (stat) and Njj = 685 ± 93 (stat). 
The shape of the Mj-^ spectrum for the sum of the 7-l-jet 
and dijet backgrounds is obtained from an independent 



control data sample by requiring Onn < 0.1 for one of 
the photon candidates, and is parameterized with an ex- 
ponential function. The resulting shape is found to be 
in excellent agreement with that derived by directly ap- 
plying the 4x4 matrix method bin-by-bin in the final 
selected sample, but has smaller statistical fluctuations, 
especially in the high My-y region. 

After subtraction of the ZDY, 7+jet and dijet back- 
ground contributions, the Mj^ spectrum is examined 
for the presence of a narrow resonance. For each as- 
sumed Mh value (between 100 and 150 GeV, in steps of 
5 GeV), the search region is defined to be {Mh — 15 GeV, 
Mh + 15 GeV), where 15 GeV corresponds to about 
five times the expected My^ resolution. The DDP back- 
ground is estimated by performing a sideband fit to the 
M^j spectrum in the 70 to 200 GeV range (this ex- 
cludes the search region) using an exponential function 
(see Fig.[Tl3). Such a parameterization has been validated 
using a next-to- leading-order for this process [l9| . 

Systematic uncertainties affecting the normalization 
and shape of the M^^ spectrum are estimated for both 
signal and backgrounds. Uncertainties affecting the ZDY 
background normalization include: integrated luminosity 
(6.1%), electron misidentification rate (14.3%) and ZDY 
cross section (3.9%). Such uncertainties are propagated, 
via the 4x4 matrix method, to the estimated normaliza- 
tion of the 7-f jet and dijet background contributions, af- 
fected in addition by the uncertainty on the Onn > 0.75 
selection efficiency for photons (2%) and jets (10%). The 
uncertainty in the shape of the 7-l-jet and dijet M^^ spec- 
trum is given by the statistics of the control data sam- 
ple used to parameterize it. The above uncertainties, 
as well as the statistical uncertainties of the sideband 
fitting method, result in systematic uncertainties in the 
normalization and shape of the DDP background contri- 
bution. Uncertainties affecting the signal normalization 
include: integrated luminosity (6.1%), acceptance due 
to the photon identification efficiency (6.8%) and PDFs 
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Mh (GeV) 


100 


110 


120 


130 


140 


150 


77 
-yj + jj 


55±7 
540±66 


17±3 

AQ^ 1/10 

319±39 


6±2 
204±25 


5±1 
OOP 1 9n 

133±16 


4±1 
89±11 


3±1 

1 0/14-00 

61±8 


total background 


1337±29 


817±26 


534±19 


374±12 


254±7 


188±5 


data 


1385 


827 


544 


357 


270 


202 


HsM signal 


1.62±0.11 


1.61±0.11 


1.51±0.10 


1.26±0.08 


0.90±0.06 


0.54±0.04 


acceptance (%) 


19.9,18.8,20.3 20.4,19.9,21.6 21.0,20.6,22.3 21.5,21.2,22.9 21.8,22.0,23.5 22.1,22.2,24.1 



TABLE I: Numbers of selected events in data, expected backgrounds, expected Hsm signal and signal acceptance (for each 
production mechanism: GF, VH, VBF), in the search region for different Mh values. The expected signal includes contributions 
from GF, VH and VBF processes, the latter two representing ~ 21 — 24% of the total signal. 



Mh 


(GeV) 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


exp. 


a X B 


88 


78 


71 


63 


58 


53 


49 


45 


41 


39 


36 


obs. 


a X B 


78 


110 


112 


76 


50 


46 


42 


53 


56 


66 


52 



TABLE IL Observed and expected 95% C.L. upper limits 
on a X B (in fb) for different Mh values. The expected 
limit is defined as the median of the distribution of limits 
in background-only pseudo-experiments. 



(1.7-2.2%) Finally, the location of the peak in the 

Mjry spectrum for signal is affected by the uncertainty in 
the relative data to MC photon energy scale (0.6%). 

Table [J shows the number of events in data, expected 
background and expected Hsm signal in six different 
search regions. The inset in Fig. [T]d illustrates the M^^ 
spectrum in the search region for Mh = 130 GeV, found 
to be in good agreement with the background predic- 
tion. The spectrum in the search region is used to 
derive upper limits on the production cross section times 
branching ratio for —> 77 (cr x S) as a function of Mh- 
The SM prediction for the ratio of the production cross 
sections for the three signal production mechanisms is 
assumed. Limits are calculated at the 95% C.L. using 
the modified frequentist approach with a Poisson log- 
likelihood ratio test statistic (2^, 21 1 . The impact of sys- 



tematic uncertainties is incorporated via convolution of 
the Poisson probability distributions for signal and back- 
ground with Gaussian distributions corresponding to the 
different sources of systematic uncertainty. The correla- 
tions in systematic uncertainties are maintained between 
signal and backgrounds. 

The resulting limits on ax B are given in Table HIl and 
displayed in Fig. [it. Although the SM Higgs boson is 
used as a reference model, the fact that the signal accep- 
tance is found to be almost independent of the produc- 
tion mechanism (see Table H]), makes the estimated limits 
applicable to other models of new physics with a narrow 
resonance decaying into 77. In the context of models of 
EWSB with enhanced B{H ~+ 77), the current search 
excludes a Hf boson with Mh < 101 GeV at 95% C.L., 
improving (slightly) upon previous results at the Teva- 
tron 22,1. While none of the other EWSB scenarios ex- 



plored can currently be excluded, the expected sensitivity 
is within less than a factor of four of the prediction for 
the model for Mh < 110 GeV, and only a factor ~ 20 
above the SM prediction for 115 < Mh < 130 GeV. As 
a result, this search contributes to the overall sensitiv- 
ity of the SM Higgs boson search at the Tevatron from 
the combination of multiple channels 23[ . Assuming the 
same integrated luminosity in all channels and a single 
Tevatron experiment, this analysis is expected to improve 
the combined upper limit on the SM Higgs production 
cross section by - 5% for 115 < Mh < 130 GeV. Fi- 
nally, this search is used to derive 95% C.L. upper hmits 
on B{H -f-f) between 14.1% and 33.9% for Mh in the 
range 100-150 GeV, in the case of models where the Higgs 
boson does not couple to the top quark. Conversely, for 
models where the GF production mode is available, this 
inclusive search allows improvement of the upper limits 
on B{H 77) to 3.4%-7.2% in the same mass range. 
These represent the most stringent limits on B{H —^ 77) 
for Mh in the range 100-150 GeV to date, significantly 
improving upon previous results from LEP and the Teva- 
tron [12]. 

In summary, we have performed an inclusive search for 
a narrow resonance with mass between 100 and 150 GeV 
decaying into 77 at the Tevatron. This channel is used 
to increase the overall sensitivit y o f the SM Higgs boson 
search program at the Tevatron [23| and allows the probe 
of new physics models predicting an enhanced rate for 
H 77. 
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